Coherent Cherenkov Radiation from Cosmic-Ray-Induced Air Showers 

K.D. de Vries,^ A.M. van den Berg,^ O. Scholten/ and K. Werner^ 

^ Kernfysisch Versneller Instituut, University of Groningen, 9747 A A, Groningen, The Netherlands 
^SUBATECH, University of Nantes - IN2P3/CNRS- EMN, Nantes, France 

(Dated: July 5, 2011) 

Very energetic cosmic rays entering the atmosphere of the Earth will create a plasma cloud moving 
with almost the speed of light. The magnetic field of the Earth induces an electric current in this 
cloud which is responsible for the emission of coherent electromagnetic radiation. We propose to 
search for a new effect: due to the index of refraction of air this radiation is coUimated in a Cherenkov 
cone. To express the difference from usual Cherenkov radiation, i.e. the emission from a fast moving 
electric charge, we call this magnetically-induced Cherenkov radiation. We indicate its signature 
and possible experimental verification. 
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INTRODUCTION 



It is a well-known fact that an emitting source moving 
at a velocity exceeding the wave-propagation velocity in 
the medium will induce a shock-wave. Prime examples 
of this phenomenon are the sonic boom emitted by a 
super-sonic airplane, a bow-wave from a moving ship, 
and Cherenkov radiation emitted by an electric charge 
moving at almost the vacuum speed of light in a medium 
with an appreciable index of refraction. In this work wc 
show that a similar effect occurs in the emission from 
a fast moving electric current. It is suggested that this 
effect manifests itself in the emission of electromagnetic 
waves from particle cascades in the atmosphere of the 
Earth initiated by ultra-high energy (UHE) cosmic rays, 
with energies in excess of 10"'^^ eV. 

An UHE cosmic ray entering the atmosphere of the 
Earth creates a cascade of particles, called an exten- 
sive air shower (EAS). In this cascade there are copious 
amounts (> 10^, depending on the initial energy of the 
cosmic ray) of electrons and positrons forming a small 
plasma cloud. This cloud, with a typical size of less than 
1 m, moves with almost the light velocity towards Earth. 
The magnetic field of the Earth, by means of the Lorentz 
force, induces a drift velocity for the leptons which is 
perpendicular to the direction of the initial cosmic ray 
and opposite for electrons and positrons. As a result 
an electric current is created in the fast moving plasma 
cloud. The strength of the induced electric current is 
roughly proportional to the number of charged particles 
in the plasma cloud as the induced drift velocity is vary- 
ing little with height. This macroscopic picture Pl was 
recently confirmed [2] to agree with a microscopic de- 
scription [3]. Even when the index of refraction of air 
would be equal to that of vacuum this varying electric 
current emits electromagnetic waves and coherent emis- 
sion occurs at a wavelength longer than the size of the 
charge cloud, i.e. for radio frequencies v < 300 MHz [3]. 
The geomagnetic emission mechanism |5] has been con- 
firmed from data [6J. In addition to the induced cur- 
rent the plasma cloud has a net charge excess which also 



radiates. The polarization direction of radiation distin- 
guishes the emission due to the charge excess and the 
geomagnetic current [7j [8] . Since charge-excess radiation 
is generally smaller in intensity we will concentrate in 
this work on geomagnetic emission. 

As is well known the propagation speed of electromag- 
netic waves is c/n where c denotes the velocity of light 
in vacuum. In this work we investigate the effect of the 
index of refraction of air, n, on the emission following 
Ref. [3]. The effects of Cherenkov radiation from EAS 
have also been addressed in Ref. [10 . In this work we 
show that for realistic values for n the Cherenkov effect 
introduces distinct features in the ground pattern of the 
emitted radiation. 



RADIO WAVE EMISSION 

As described, a cosmic ray entering the atmosphere in- 
duces an EAS, a plasma cloud moving at almost the light 
velocity, where the magnetic field of the Earth induces 
a net electric current in the plasma [T]. From classical 
electrodynamics |llj we obtain the Lienard-Wiechert po- 
tentials for a point source following a trajectory ^(t') and 
an observer at rest at {t,x), 






{t,x-^{t')) = ^^ 



(1) 



where f denotes the retarded time corresponding to the 
time the signal was emitted from the moving charge dis- 
tribution and V is the retarded distance [S] . In the point- 
like approximation (denoted by pl) where the size of the 
plasma cloud is infinitesimal the four-current is. 



rp^it',x) = j'^{t')s'{x-m), 



(2) 



where f (t') — — ci'e^ and JpL carries a longitudinal com- 
ponent due to the net charge excess (which we will ignore) 
and a transverse component due to the opposite drift 
of electrons and positrons induced by the Earth mag- 
netic field. The transverse component is proportional 
to the number of charged particles at a given height. 



J-'-(t') ex Nf,{z) where z = —ct' is the distance to the 
Earth's surface as measured along the shower path. 
The electromagnetic fields at the observer are given as, 



E' = c(aM° - d°A') , 



(3) 



with i = x,y, z, an observer located at x — {x, y, z — 0), 
and where at t = t' = the shower hits the Earth. The 
distance d = ^Jx^ + y^ is equal to the transverse dis- 
tance to the shower axis. The velocity of the charge 
cloud is c/3 where we set /3 = 1 in the following. For 
the special case n = 1, the retarded distance T) — 
^{—fictY + (1 — "n?-^'^)(P reduces to I? = /3ci which is 
finite. 

The retarded time is defined by the light-cone condi- 
tion c{t—t') = L{x, £,{t')) , where the distance L = L{x, ^) 
is the optical path length between the source located at 
^ and the observer at x. In reality the index of refrac- 
tion of air depends on density and thus height, n = n(z), 
hence light will follow curved trajectories where L is the 
integral / ds along the light curve from ^ to x. 
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FIG. 1: The emission height z = —ct' as function of the 
observer time t for three different values of the index of re- 
fraction. The dashed line gives the shower-profile as function 
of z for a E — 5 X 10^^ eV proton-induced shower. 

In Fig. [l] the emission height, z — —ct', is plotted as 
a function of the observer time t for an observer 100 m 
from the shower axis for three choices for the index of 
refraction n. For the case n = 1 the plot of the retarded 
time (red drawn curve) shows that the retarded time is a 
single valued function and that the earliest signals come 
from the top of the shower. For an index of refraction de- 
viating from unity (n = 1.0003, n — n{z)) the function is 
composed out of several branches (dashed-magenta and 
dotted-blue curves in Fig. fl]) limited by the critical times 
t = tc, where the derivative dt' /dt becomes infinite. In [9] 
it was already shown that a singularity in this derivative 
is related to a singularity in our vector potential since 
V — R/{dt' /dt). This singularity (the branch point) 



is well known and corresponds to Cherenkov emission. 
Since the critical time is the time the Cherenkov radia- 
tion is seen by the observer, it is henceforth called the 
Cherenkov time. In the case of a constant and finite re- 
fractivity N — n — 1 the Cherenkov time tc is given as. 



etc = d \J(v?- - 1) , 



(4) 



with corresponding retarded time —t!^ — tdir? — 1). The 
retarded time can readily be converted into an emission 
height Zc = —ct'c and it can be seen that the expres- 
sions are consistent with the expression of the angle for 
Cherenkov emission, cos 0c = ^jn. For the general case 
it has to be calculated numerically. 

To calculate the realistic pulse form for n 7^ 1 it is 
essential to include the fact that the emitting plasma 
cloud has a finite size, 

Al (t, f ) = Jd^rjdh p{h, r) A^'p^ (i, f - () (5) 

where r — rxe^ ~\- rity is the relative transverse coordi- 
nate to the shower axis and the source is at the position 
£_{t',h,ri) = [—ct' + h)ez + f. The density profile is 
parameterized as [TH [13] 

w(r, h) = 27rrp(/i, r) = N{r + ro)"^-^/i e-^'V^ (g) 

at a fixed shower time t' and where h is the longitudinal 
distance from the shower front that, by definition, moves 
with the vacuum speed of light c. The normalization 
constant, N, is chosen such that Jdrjdhw{r,h) — 1. 
Positive values of h mean a position behind the shower 
front, hence w is zero for values of /i < 0. The parameters 
are taken as rg = 80 m, L = 0.5 m following the results of 
simulations using the cascade mode of the CONEX-MC- 
GEO shower simulation and analysis package [3] . The re- 
sults of these simulations indicate that the pancake thick- 
ness parameter L has to be considerably smaller than 
the value of L = 3.9 m used in previous calculations [7]. 
Retaining the terms pertinent to geomagnetic radiation 
and using partial integration to have the derivatives in 
Eq. (|5| operate on the density distribution exclusively, 
the electromagnetic field can be expressed as 



E' =- 



Mqc 

in 



d^r 



dh 



1 

\V\ 



9w „ 



dJ 



PL 



dt' 



(7) 



The complication in Eq. ([7]) is now reduced to the inte- 
gration of inverse square-root divergencies over smoothly 
varying functions giving finite results. We will argue later 
that the first term is important for Cherenkov radiation 
while the second is dominant for n = 1. 

To get a better appreciation of the structure of the 
pulse and the influence of the index of refraction we will 
discuss in some detail the results for an observer at a dis- 
tance d = 100 m from a shower with E = 5x 10^^ eV at an 
angle 6 = 27° . To simplify the discussion we will -for the 
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FIG. 2: Pulse due to geomagnetic radiation at d = 100 m 
from the core for different values of the index of refraction as 
discussed in the text, n = 1, n = 1.0003 fixed, and n = n{z) 
realistic. In this calculation the lateral extent of the shower 
is ignored. 



time being- ignore the radial extent of the charge cloud. 
Essential for the structure of the pulse is the longitudinal 
shower profile Ne{z), as shown by the thinly-dotted black 
curve in Fig. IT] 

For the case that n = 1 the denominator I? is a 
smoothly changing function over the integration regime 
and the dominant contribution is derived from the last 
term in Eq. ([T]) . Since the derivative of the shower pro- 
file reaches a maximum at z = 6 km it can be seen from 
Fig. [T] that the peak of the pulse occurs at i « 3 ns which 
agrees with the result in Fig. [2J The integral of the first 
term, being a derivative, almost vanishes. 

For the cases that n ^ 1 the reasoning is rather differ- 
ent. The denominator V vanishes at the Cherenkov time, 
indicated by the vertical arrows in Fig. Ill and thus l/V 
varies strongly as function of h. The contribution from 
the first term in the integral Eq. ^ is large (in contrast 
to the case n = 1) and results in an enhanced contribu- 
tion from the corresponding emission height. The pulse 
height will thus be proportional to Nf,{zc) while the peak 
is observed a little after t^ in agreement with the re- 
sults shown in Fig. |2] Since J^{t'^) oc Nf,{zc) is large for 
d = 100 m this results in a strong pulse. The remaining 
contributions in the integral are of secondary importance 
in this case. 

In reality the refractivity is equal to iV(0) = 3 x 10^'' 
at ground level and decreases exponentially with height, 
N{z) = n{z) — 1. The values of the retarded time thus 
lie in between those obtained with n = 1 and n — 1.0003 
and are shown as the dotted-blue curve in Fig. [T] Also 
for this case there is a clear Cherenkov time which is 
slightly smaller than that for n = 1.0003 resulting is a 
very similar pulse as for the case of a fixed-finite n. 

The main effect of including the radial extent of the 
shower is to smooth the time structure of the pulse and 
thus to wash-out some of the effects of n 7^ 1. The third 



panel of Fig. [3] gives the complete geomagnetic result and 
should be compared with Fig. [2] The differences between 
the three different choices for the index of refraction have 
diminished, instead of being three times as large the pulse 
height is increased by a factor two only. 
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FIG. 3: Pulse due to geomagnetic radiation at at several dis- 
tances from the core for different values of the index of refrac- 
tion as discussed in the text (n = 1, n = 1.0003 fixed, and 
n — n{z) realistic), including the lateral extent of the shower. 

The Cherenkov effect is strongly distance dependent as 
shown in Fig. [31 As expressed in the discussion follow- 
ing Eq. Q emission at the Cherenkov angle cos^c = l/n 
relates a distance d from the shower core to a particular 
emission height. For an observer close to the shower, 
d = 20 m, this Cherenkov height lies well below the 
shower maximum where Ne{z) has fallen considerably 
and the intensity of the emitted radiation is low. It 
increases with increasing distance to reach a maximum 
around d — 100 m where the Cherenkov point lies at the 
height of the shower maximum. At even larger distances 
the Cherenkov point lies above the shower maximum and 
the amplitude decreases rapidly. The position of the 
Cherenkov maximum thus reflects the shower maximum 
and is thus sensitive for the composition of the cosmic 
ray. For n = 1 the main strength of the pulse is emitted 
at the height where the derivative of the shower profile 
is large. Since the retarded distance, |I?|, correspond- 
ing to this height in the shower evolution increases for 
observers further away from the shower core, the ampli- 
tude of the pulse is a monotonically decreasing function 
with distance. 

Calculations confirm that for charge-excess radiation 



(not reported here) the effect of tlie index of refraction is 
very similar as for geomagnetic radiation where there are 
subtle differences due to the somewhat different weight- 
ing over shower height. For the chosen geometry the ge- 
omagnetic effect is however dominant. This shows that 
the Cherenkov effect applies equally to radiation from a 
moving charge, to which it is usually applied, as to that 
from a moving electric current, which is at the focus here. 



EXPERIMENTAL VERIFICATION 



Some hints of this effect may have been seen in re- 
cent results from the LOPES [14] collaboration show- 
ing that for certain events the pulse height follows the 
trend shown by the dotted-blue line in Fig. |4] More de- 
tailed measurements are necessary where full attention is 
given to polarization observables which are crucial to dis- 
tinguish geomagnetic and charge-excess radiation. Such 
measurements are planned for new and renovated set- 
ups at LOPES, CODALEMA [TS], and recently also new 
set-ups at the Pierre Auger Observatory (MAXIMA [TB] , 
AERA [T7]), and LOFAR [TH]. 



The polarization of the radio signal, which distin- 
guishes geomagnetic and charge-excess radiation, is not 
affected by the Cherenkov effect. The principal signa- 
ture of the Cherenkov effect is that the pulse at a cer- 
tain distance, 100 m for the present example, is consid- 
erably larger than it would have been for the case that 
n = 1. This feature is especially clear from Fig. [4] where 
the height of the pulse is plotted at various distances from 
the core. At short distances the pulse height diverges for 
the case of n = 1 where it should be noted that this di- 
vergence is strongly dependent on the pancake thickness 
parameter L, smaller values give a stronger divergence at 
d = 0. Apart from an overall scaling factor the picture 
is not affected by L ior d > 50 m. For the realistic case 
where n ^ 1 a marked deviation is predicted with a clear 
enhancement in peak intensity at distances ranging from 
50 till about 100 m from the core. This peak results from 
the fact that the Cherenkov point lies close to the shower 
maximum for d =100 m. 
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FIG. 4; Pulse height as function of distance to the shower core 
for the three choices for the index of refraction as discussed 
in the text {n — 1, n — 1.0003 fixed, and n — n{z) realistic). 

The feature shown in Fig. |4J an increase of the ampli- 
tude of the geomagnetic radiation with distance for n ^ 1 
instead of the monotonic decrease predicted for n = 1, 
remains for showers at a large angle with the vertical. 
For these inclined showers the shower maximum will oc- 
cur at a larger values of z, the distance to the point of 
impact on Earth, and the maximum in the intensity will 
thus be observed at larger distances from the shower core, 
consistent with the angle for Cherenkov emission. 
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